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I. INTRODUCTION 



The relationship between the ability to split off the centre of mass motion, the idea of 
a "guiding centre" and its connection with some form of generalised Galilei, or Newton- 
Hooke type " kinematic symmetry" [TJ |2] has been the subject of a number of recent papers 
[3H5]. In the presence of magnetic fields this is the subject of Kohn's theorem [6j and its 
variants. The use of the guiding centre approximation in plasma physics is well known. Less 
well explored is the application of these ideas to gravitational physics. It is true that the 
idea of a guiding centre is well established in galactic dynamics [7] , but its connection with 
kinematic symmetries does not appear to have been explored before. The purpose of the 
present paper is to fill that gap in the literature. 

The oldest example of what we have in mind are Hill's equations for the Earth- Moon-Sun 
system [HI [9] . However with the development of our understanding of the structure of the 
galaxy, it was realized that similar equations hold for the motion of stars around the Milky 
Way [7J [TUHT2]. Understanding many electron atoms in the Old Quantum Theory leads 
to the same equations and its failure to deal with the Helium was the notorious stumbling 
block which led to the development of modern Quantum Mechanics. In more recent times 
there has been a revival of interest in semi-classical models of many electron atoms [13] and 
to muonic atoms [14j . 

The plan of the paper is as follows. In section ITT] we introduce and derive Hill's equations. 



In section III we analyze their symmetry group and its relation to the centre of mass motion. 



and in particular (in the planar case) show that it is five-dimensional. In section IV we obtain 
the Lie algebra using its vector field generators acting the Newton-Cartan spacetime. In 
section [V] we pass to a Hamiltonian treatment and show that the Poisson algebra of moment 



maps is an extension by two central elements. In section VI we provide the Eisenhart-Duval 
[ToTHS} |2"T] lift of the system to a 3 + 1 dimensional metric with Lorentz signature which is not 



conformally flat, as we show explicitly. In section VIII we give an alternative interpretation 



of the Hill system in terms of a Landau problem in an anisotropic oscillator and in section IX 



we use this representation to give a "Chiral Decomposition" using the methods of [221 123] . 
Section [X] describes some possible variants and extensions of our results and the last section 
is a short conclusion. 
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II. HILL'S EQUATIONS 



As a model for the Earth- Moon-Sun system [HI E] , or for a cluster of stars moving around 
the galaxy in an approximately circular orbit [7J [TUHT^] . one has the following equations 

i (x - 2uv - 3u 2 x ] - V Gm - m ^ ~ x ") 



m„ 



a {y a + 2ux a ) = ^ u-gja ' (1) 

b=£a 

, : \- Gm a m b (z b - z a ) 

a (z a + u 2 z a ) = ^ 
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These equations are valid in a suitable rotating coordinate system. For the Earth- Moon-Sun 
system, for example, X\ can be the position of the Earth and x 2 can be that of the Moon. 
The motion of the Sun is neglected, and the remnant of its influence on the Earth-Moon 
pair is represented, in the first-order approximation, by the repulsive anisotropic harmonic 
term [36J in the first equation, where 

is the angular velocity of a circular Keplerian orbit lying in the z = plane and having 
radius R. The linear-in-velocity terms correspond to the Coriolis force induced in a rotating 
coordinate system. The right hand sides represent the gravitational interactions between 
the Earth and the Moon. 

These equations are obtained as follows (221 EZ]- Let r,9',z be cylindrical coordinates 
centered on the Sun or on the galactic center which has a mass M so large compared with 
those of the other moving masses that it may be assumed to remain at rest in an inertial 
coordinate system. Let x, y, z be coordinates with respect to a rotating coordinate system 
whose origin lies on the a Keplerian orbit with r = R, 8' = cut, z = 0. The x axis is taken 
to be radial so that r = R + x , and the y axis is taken to be tangential to the orbit. The 
forces acting on each particle whose comoving coordinates are x a consist of their mutual 
gravitational attractions and the attraction due to the gravitational potential produced by 
the Galaxy 

u=-y , GMma (3) 

To quadratic accuracy in x a ,y a , z a 

^ GMm a / x a 1 y\ + z 2 a x\ \ 

^ R \ R 2 R 2 + RV' {) 



which implies that the force is to linear accuracy 

- VC/ = ^2m a u 2 (-R + 2x a ,-y a ,-z a ) . (5) 

a 

Substitution in Newton's equations of motion now gives Hill's equations Q. 

Strictly speaking the equation originally considered by Hill was a special planar case for 
two bodies in which the position of the Earth xi was assumed to be at rest X\ = and the 
position of the Moon X2 thus to satisfy 



Gm\X2 



y 2 + 2cux 2 



Gm x y 2 



We omit henceforth the z variables and work in the plane. We mention however that the 
more general case where motion in the z direction is allowed has important applications, 
either to the Earth- Moon-Sun system, see [9] , or in semi-classical treatments of the helium 
atom [13] or muonic atoms [H]. 



III. SYMMETRIES AND CENTER OF MASS MOTION 

In addition to the discrete symmetries of parity and time reversal, 

x a (t) ->■ -x a (t) and (x a (t), y a (t)) (x a (-t), -y a (-t)) : (6) 

Hill's equations admit a continuous four-parameter family of abelian symmetries, since they 
are invariant under "translations and boosts" 123] , 

x a ->■ x a + a(t) . (7) 

Inserting into the Hill equations and putting (with some abuse of notations) a = (x, y), 
allows us to infer that to be a symmetry requires 

x — 2uii — 3u 2 x = 

(8) 

y + 2cox = 

The simplest way to solve these equations is to derive the upper equation w.r.t. time and 
then use the lower equation to eliminate y to yield an oscillator equation for x, cPx/dt 2 = 
—uj 2 x. Thus 

, , A , B 
x[t) = — smut cosut + xq. 
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FIG. 1: Trajectory of the center of mass in the Hill problem. The straight horizontal line in 
the middle indicates the trajectory of the guiding center about which the center of mass performs 
"flattened elliptic motion". 

Putting x(t) into the second equation and integrating provides us with y(t); testing the pair 
x(t),y(t) on our original system fixes the integration constants to yield 



Then 



A . B 

xyt) = —smut coswt + xq 

UJ UJ 

. A B 3 

y(t) = 2 — coswt + 2 — sinwt cotxo + yo. 

uj uj 2 



lo 2 (x-x ) 2 u 2 (y - y + lux ty 



A 2 + B< 



+ 



A(A 2 + B 2 



(9) 



(10) 



shows that the trajectories are ellipses centered at (x ,y — ^cjx t) with major axes lying 
along the y direction. The ratio of the semi-major to the semi-minor axis is 2 :: 1, and the 
centers drift along the y direction with constant speed |cjxo in the direction of its major 
axis, see Fig. [Tj 

At this point, it is legitimate to wonder what is the interest of studying the properties of 
(pj), which only describe some property of the system but not the physical system itself. A 
justification comes from observing that, as a consequence of the linearity of the l.h.s. of 
and because the gravitational forces on r.h.s. of these equations satisfy Newton's third law, 
the center of mass, 

X | _ Y;a m a X a 

y 



X 



(11) 



(with another abuse of notation) satisfies exactly the same equations ([8]). The latter de- 
scribes therefore more than a "property". 

We emphasize that our statement relies on the equality of the inertial and passive grav- 
itational mass, m a , for objects with significant self-gravitation. In other words, both the 
center-of-mass decomposition and Galilean symmetry depend on the so-called Strong Equiv- 
alence Principle jH [18]. It has been verified experimentally by lunar laser ranging to very 



high accuracy using the Nordtvedt effect for the Sun-Earth-Moon system. If 

^passive -, Eq 

= 1- VN - (12) 

^inertial ^ ^inertial 

where Eq is the Gravitational self-energy, then 

\rj N \ < (4.4 ±4.5) x 10~ 4 . (13) 



Since 2 E ° — = 4.6 x 10 10 for the Earth and 2.1 x 10 10 for the Moon, the strong equiv- 

C Tiinertial 

alence principle is satisfied to better than one part in 10 13 [28] . 

Below we focus our attention at the symmetry alias center-of mass equation Q. 
The general solution ^ is composed of two particular cases. 

• Let us choose first xq = yo = 0; then the trajectory is an ellipse centered at the origin, 
and oriented along the y direction, 



X + (t) 



xi(t) 
xlit) 



/ A . B \ 
I — sin cot cos cot > 

CO CO 

A B 
\ 2— coswt + 2 — smcot I 



(14) 




CO CO 

Putting A = B = provides us instead with 

*-(*)=! i = f 3 Xo i- ( i5 ) 

-%Lotx + y 

The particular form of this solution comes from a delicate balance between the harmonic 
and the inertial forces which precisely cancel, 

3w 2 5 H X! + 2coe ij Xi = 0, (16) 

so that the particle drifts perpendicularly to the harmonic field with constant velocity. 
Anticipating what comes below in Sect. VIII[ we call it a Hall motion. 



As it will be explained in Sect. IX, the first of these particular solutions, namely X , 



describes the guiding center, and the second, X + , describes the relative motion around it. 
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IV. VECTOR FIELDS AND ALGEBRA 



The planar symmetry group is generated by the spacetime vector fields 

K\ = — ( sin utd x + 2 cos utdy) 
Kl = — ( — cos utd x + 2 sin u>td v ) 

Kl = d X ~ \utdy ( 17 ) 
Kl = dy 

H = d t . 

Here the vector fields K±, i — 1,2 generate the infinitesimal time-dependent symmetries 
([7]), and H represents infinitesimal time translations. The non-trivial brackets are 

[H,K\] = -uK%, 

[H,Kl] = +uKl, (18) 
[H,Kl] = -\uKl. 

Normalizing the total mass to unity, Y2a m a = 1' the Lagrangian for the center of mass 

is 

1 3 

L = 2 ( ±2 + y 2 ) ~ ~ x y) + 2 w2a;2 ' ( 19 ) 

The mechanical momenta p x = x and p y = y do not Poisson-commute, {p x ,p y } = 2w . The 
Hamiltonian is 

tf=i(g+p})-|*V. (20) 
It is Liouville integrable since H, and the dual momentum p y = y + 2oox mutually commute 

V. MOMENT MAPS FOR HILL'S EQUATIONS 

Following [3J, we find the conserved quantities 

1 1 / \ 
ki_ — — ( sin cut + 2p„ coswt + 3wa; coswt), 
u; 

k+ = — ( — p x cos cut + 2p v sin cot + 3ux sin cut) , 
3 

«i = {Px - - 2 Wt (Py + 2wx )' 
k 2 = + 2cux . 



(21) 



Recovering the generating vectorfields in (17) as 



- K = {k, x}d x + {k, y}d y 



can be viewed as a consistency check. 



Note that the Poisson algebra does not coincide with the Lie algebra (18) since k\ and 
k 2 , and k\ and k 2 _, do not Poisson commute but their brackets give rather two central 
extensions, 

1 (22) 
{kI,k 2 _} = --U. 



The other bracket relations (18) are unchanged. Thus the conserved quantities realize two 
commuting copies of Heisenberg algebras. Evaluating the moment maps on the solutions 
@ gives 



i _ B 
A 



^ (23) 



«i = —^yo, 

2 1 

which shows that they are indeed constants of the motion. Hence n 2 _ = \ojxq commutes 
with the Hamiltonian H, and is related to the [conserved] x coordinate of the center of the 
revolving ellipse. In terms of the conserved quantities, the Hamiltonian reads 

H = y (K) 2 + (4) 2 ) " \{#- ) 2 = \(A 2 + B 2 ) - \^xl . (24) 



One thus has 



{H,k\} = —uk 2 + 



{H,k 2 + } = +uk\ 

{H,Kl} = -\u>#_ 

{H,k 2 _} = 0. 



(25) 



The k, although explicitly time-dependent, are however conserved, — = — + {«, if} = 

(JjI/ yJv 
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VI. EISENHART-DUVAL LIFT 



Following the procedure described in [3J [TBI HE] we lift the Hamiltonian to that of a 
massless particle in 3+1 spacetime dimensions. The calculation is straightforward and we 
just give the result. The 4-metric is given by 

ds 2 = dx 2 + dy 2 + 2dt(dv + u(xdy-ydx))+3u 2 x 2 dt 2 (26) 

where v is a new, "vertical" coordinate [16]. 



As pointed out in Ref. [18] . the 2cu-terms in (26) admit a two-fold interpretation. In 
the present context here, they can be viewed as representing inertial forces in our rotating 



coordinate system. In Section VIII below they will be interpreted as an external magnetic 
field. Its null-geodesics project onto "ordinary" spacetime according to the center-of-mass 
[alias symmetry] equations of motion (j8j). 



(26) is a Ricci-flat 3+1 dimensional Lorentzian metric with covariantly constant null 
Killing vector field 

£ = d v , (27) 

i.e., a "Bargmann space" [16]. 

The "Bargmann" framework is particularly convenient for describing the symmetries. 



Our symmetry transformations lift indeed to the "Bargmann" metric (26) as isometries. 
Let us assume that ([7j satisfies the symmetry condition ([§]). Completing ^ with 

v — > v — i (a ■ a + 2a • x + 2cu a x x) , (28) 



a tedious calculation shows that the "Bargmann" metric ([26]) is left invariant. 
Working infinitesimally, the vector fields (fT7|) lift as 



K\ = — ( sin cut d x + 2 cos cot d y ) + (x cos cut + ysmcot)d v , 
cu 

K\ = — ( — cos cut d x + 2 sin cut d y ) + [x smut — y cos ut)d v , 
cu 

~ 3 3 1 

K l _ = d x - -cutdy + ( - -cu 2 xt + -cuy)d v , 

K 2 _ = d y + ouxd v , 

whose Lie brackets are found to be 

{K\,Kl} = -U 

cu 
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(29) 



(30) 



which are, up to sign those, (22), satisfied by the associated conserved quantities [37] . 

The lifted symmetries realize hence not original Lie algebra structure, (18), but rather 
their central extension with £, the generator of vertical translations, as central element. 



VII. BARGMANN SPACES WITH NEWTON-HOOKE SYMMETRY 

The origin of Newton-Hooke symmetry has been understood a long time ago [IB] : the 
Bargmann space of an isotropic harmonic oscillator with time-dependent spring constant 
k(t) is 

dx 2 + dy 2 + Idtdv - k(t)(x 2 + y 2 )dt 2 , (31) 

and the massless dynamics "upstairs" projects on the oscillator dynamics "downstairs". 
Newton-Hooke symmetry is represented by the isometries of this metric, and is indeed a 
subgroup of £ = 5,, -preserving conformal transformations — the latter forming the (centrally 
extended) Schrodinger group. 



The metric (31) is, furthermore, Bargmann-conformally flat i.e. can be mapped confor- 
mally onto 4d Minkowski space by a ^-preserving transformation. 

Bargmann-conformally related metrics share the same symmetries; a conformally flat 
Bargmann metric admits therefore the same [namely Schrodinger] symmetry as a free par- 
ticle. 

Now we describe, following Refs. [TT j IT9H2T] . these Schrodinger-conformally flat spaces. 
lnD = d + 2>3 dimensions, conformal flatness is guaranteed by the vanishing of the 



conformal Weyl tensor 



' & D"] i 2 xli* X"\ 



= «"„ - — 2 ^iT< + (C _ 1)(J _ 2) VI, «■ (32) 

Now R^vpa^ = for a Bargmann space, implying some extra conditions on the curvature. 
Inserting the identity = into C^ pa = 0, using the identity ^R% = {R v a = R^ a ), 

we find 

o = - [z p k - z*r;] + - ^1 • 

Contracting again with £ CT and using that £ is null, we end up with R^, p ^ u = 0. Hence the 
scalar curvature vanishes, R = 0. Then the previous equation yields = an d thus 

R v a = ^rf for some vector field rj. Using the symmetry of the Ricci tensor, R\p V ] = 0, we 
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find that 77 = for some function g. We finally get the consistency relation 



R 



(33) 



The Bianchi identities (V^-R^ = since R = 0) yield ^d^g = 0, i.e. g is a function on 
spacetime. The conformal Schrodinger-Weyl tensor is hence of the form 



pa 



pa 



D 



(34) 



It is noteworthy that Eq. (33) is the Newton-Cartan field equation with g/(4irG) as matter 



density. Eq. (33) also implies that the transverse Ricci tensor of a Schrodinger-conformal 



flat Bargmann metric necessarily vanishes, Rij = for each t. 

Further results are only worked out for total Bargmann dimension D = 4. Since the 
transverse space is d — 2-dimensional, R^ = implies that the latter is (locally) flat and 
we can choose gij = gij{t). Then a change of coordinates (x,t,v) — > (G(t) x,t,v) where 
G = (Gij) is the square- root matrix 5 a bG^G b j = gij, casts our Bargmann metric into the 
form 

dx 2 + dy 2 + 2dt [dv + A ■ dx] - 2Udt 2 . (35) 
Now we turn to determining all such conformally flat 4-metrics. The non-zero components 



of the Weyl tensor of (35) are 



C 



xyxt 



-G 



ytts 



G 



xyyt 



+G. 



Xtts 



fxB, 
\d y B, 



C xtxt = -\ [d t (d y A y - d x A x ) - A x d y B, ] + \ [d 2 x - d 2 ] U, 



G 



ytyt 



1 [dt(d y A y - d x A x ) - A y d x B] - \ [dl - d 2 ] U, 



C 



xtyt 



[d t (d x Ay + d y A x ) + 2d x d y U] - \{A x d x - A y d y )B. 



Then Schrodinger-conformal flatness requires 



\eijB{t)xi + aj 




V x a = 0, 



d t a = 0, 



\C(t)r 2 + F(t) -x + K(t). 



(36) 



Note, en passant, that Eq. (33) automatically holds in this case, because 

Rxt — 2C x tts — , R 



yt 



^Cytts — 0. 



(37) 
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The only non- vanishing component of the Ricci tensor is 



Rtt = -<9*(V • A) - l -B 2 -AU= ~\B{tf - 2C{t). (38) 



The metric (36) describes a uniform magnetic field B(t), an [attractive or repulsive, 
C(t) = ±u 2 (t)] isotropic oscillator and a uniform force field F(t) in the plane. All 
fields may depend arbitrarily on time. It also includes a curl free vector potential a(x) 
that can be gauged away if the transverse space is simply connected: a, = dif and the 
coordinate transformation (t,x,v) — > (t,x,v + /) results in the 'gauge' transformation 
Ai —> Ai — dif = —\B e,ijxK If, however, space is not simply connected, we can also in- 
clude an external Aharonov-Bohm-type vector potential, explaining the o(2, 1) conformal 
symmetry of a magnetic vortex [29] . 

Our "one-sided" anisotropic oscillator here does not qualify therefore: its Weyl tensor 
does not vanish due to the anisotropy: 

[d 2 x -d 2 y ]U^0 for U(x,t) = -^u 2 x 2 . (39) 



The 4-metric (26) can not be mapped conformally to empty Minkowski space. 

How can it have the "Newton-Hook-type" symmetry, then ? There is no contradiction, 
though. Let us stress that we did not find here full Newton-Hooke symmetry, only its time 



dependent translational part : rotational symmetry is plainly broken for the metric (26). 
The latter does not come therefore by "importing" from the free case. 



VIII. RELATION TO THE LANDAU PROBLEM 

Now we point out that the center-of-mass Hill system (j8l) can also be be viewed as a 
charged anisotropic harmonic oscillator in a uniform magnetic field described by the planar 
Hamiltonian system 

{*V}=0, {,•'./,'} {j?,pl}=eBe*, H = ^ + ^(x 1 ) 2 + ^(x 2 ) 2 , (40) 

where we still scaled the total mass to unity. Comparing the equations of motion £ = {£, H} 
implying x l — eBe^x^ + kix % = [no sum on % in the last term], with ([8]) shows that the 
Hill system can indeed be viewed as a repulsive anisotropic oscillator in a uniform magnetic 
background with 

kx = k— -3u 2 , k 2 = 0, eB = 2u. (41) 
13 



The identity of the two systems relies on the equivalence of the inertial Coriolis force in a 
rotating frame with the Lorentz force due to an external magnetic field [TTl ITS] . Note also 



that the condition (16) is then in fact the Hall law, 



■ E j 

Xi= e v— (42) 

with the identifications eE l = 3u 2 5 u X}_ and eB = 2u>. 

Let us also stress that the possibility of decomposing the magnetic field plus oscillator sys- 
tem into center-of-mass plus relative motion depends on Galilean [U |30] (more precisely, on 
Newton- Hooke El [3T] ) symmetry, which requires in turn the Kohn condition charge/mass 
= constant to hold. Furthermore, as "charge" equals "mass" here, Kohn's condition is 
automatically satisfied, providing us with the required Galilean symmetry. 



In Section IX, the system will be further analyzed by decomposing (\8ti into chiral com 



ponents along the lines of [22] as adapted to the Landau problem 



IX. CHIRAL DECOMPOSITION OF THE HILL SYSTEM 

The problem can further be analyzed by decomposing our magnetic field + anisotropic 
oscillator \k\ — k, k% — 0] system into chiral components, generalizing the trick of Ref. 
[22] [23] . Define indeed the two planar vectors X± = (X±) as [38] 

p 1 = a + Xl + a_X!, p 2 = -p+Xl- P-Xl, X = X + + X_, (43) 

where a± and (3± are suitable coefficients to be found. The symplectic form 

eB .. 

Q = dp 1 A dx l + — e l W A dx 3 , (44) 



whose associated Poisson bracket is (40), is written as 

Q = (- a+ -P+ + eB) dX\ A dX\ + (-«_ - P- + eB) dX l _ A dX 2 

+ | (-a- -p+ + eB) dX\ A dX 2 + (a+ + p- - eB) dX\ A dX v _ j . 
The symplectic form splits into two uncoupled ones when 

a_ + P+ = eB, a + + (3- = eB. (45) 
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The Hamiltonian becomes in turn 



1 k 
H = - {a 2 + X 2 + X 2 + + P%X\X\ + a 2 _ X 2 _X 2 _ + 01X^X1) + - (X\X\ + X l _X l _) 

+ | {a + a_ ) X\X 2 _ + (/?+/?_ + k) X\ X x _ } , 
which splits into H = H + + H when 

= 0, P+/3- + k = 0. 

Since our formulae are symmetric in a + and a_, we can choose a_ = to find 

k „ fc 



a_ = 0, 



/?+ = eB, (3. 



eB' 



With such a choice we will have decomposed our system as 

Q, = -( e B + — )dX] AdX 2 , + (eB + — )dX 1 A dX 2 
v eB' + + v eB' ~ 



efi(e5 + A )x i X i +(ei?+ A )2xW 



+ 



2e5 



eB 



X l X l 



(46) 



(47) 



(48) 



(49) 



#4 



Note that f2+ and f2_ have opposite signs. 



Returning to the Hill problem, inserting the matching coefficients (41) into (48)-(49) 
yields 



and hence 



p 1 = 



p 



-2ooX] - -ooX 1 . 
+ 2 - 



X = X+ + X 



I - -u dX\ A dX\ | + | ho dX x A dX 2 } , 



(50) 



(51) 



(52) 



The Poisson brackets associated to this symplectic form show that both sets of coordinates 
X\ and X l _ are non-commuting, 



{xl,x 2 + } = l {xl,x 2 } = {xl,xl} = o, {x^x 2 _} = - 2 - ) 

oo oo 



(53) 



and provide us with the separated equations of motion 



X 1 



X 2 
15 



-2uX\, 



-coX 1 . 



(54) 



whose solution allows us to recover (14)-(15) once again. The general solution d9|) is the sum 



of the chiral components, X(t) = X + (t) + X_(t). 



Here the simple X dynamics is that of Hall motion with constant velocity drift (15) 



of that of the guiding center, while the X + system, whose trajectories are those flattened 



ellipses in ( 14 ) describes the anisotropic oscillations of the center of mass about the guiding 
center. 

Having decomposed the center-of-mass alias time-dependent translation-symmetry equa- 
tion into chiral components, the Newton-Hook symmetry plainly follows from those of our 



chiral solutions. For the separated equations (|54|) the inital conditions, 
X+(0) 



1 -B/u^ 


H 


v 2A/u J 





\ 



X_(0) 



Xq 



yyoj 



Xl(t) cos ut - \X\{t) sin ut 
2X\ {t) sin ut + X\ (t) cos ut J 

\ 



(55) 



Xi{t) 



{x*(t) + -utX 1 _(t),) 



(56) 



are plainly constants of the motion. They are in fact proportional to those in Eqs. (21) and 



(23), 



X+(0) 



X. SOME VARIANTS 



-4 

-2k1 



AL(0) 




(57) 



The ideas of the present paper may be generalized in various directions and and even 
applied to areas beyond the realm of classical gravity to quantum semi-classical treatments 
of quantum systems. In this section we we briefly outline some examples. 



A. Anisotropy 

The application of Hill's equations to galactic clusters was first suggested by Bok pU] and 
by Mineur [TT] and developed by Chandrasekhar [12] , Chandrasekhar did not assume that 
the gravitational field of the galaxy was just a simple monopole. As a result he obtained 
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equations of a more general form: 

Gm a m b (x b - x a ) 



m a (x a - 2uy a - 3ufx a j = ^ 



13 



x a -x b 

of a 

,.. . . \ Gm a m b (y b -y a ) 
m a (y a + 2ux a ) = ^ ■ _ [3 (58) 

b=£a 

(.. 2 \ Gm a m b (z b - z a ) 

of a 



1 3 



where the z coordinates were restored. One still has the abelian symmetry ^ but ^ 
becomes 

x — 2uy — 3cvfx = 

y + 2cox =0 (59) 



z + <jj\z = 



and d9|) is replaced by 



A B 
x = — sin fit — — cos fit + xq, 

to to 3c<j? 

y = 2A— cos fit + 2B— sin fit - ——x t + y , 
ir \V 2lo 

z = C cos oo 3 t + D sin co 3 t , (60) 



where fi = a/4w 2 — 3cof is called the epicyclic frequency and often denoted by k. The ellipses 
now have ratio of major to minor axis equal to and move with speed ^7^0- 
The symmetry is generated by the vector fields 

to 

K\ = sin fit d x + —2 cos fit d y 
K 2 = 2 — sin fit d y — cos fit d x 

K 3 = d x - ^-tdy 
2oo y 

K, = dy (61) 

K 5 = cos u 3 td z 
K 6 = smu 3 td z 
H = d t 
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whose non-trivial brackets read 

[tf.tfi] = -QK 2 

[H,K 2 ] = +SIK X 

[H,K 3 ] = (62) 

[H,K B ] = -u 3 K 6 

[H,K 6 ] = u 3 K 5 . 

By rescaling the generators, the sub-algebra they span may be seen to be independent of 
the parameters uj and lo\. 

B. Electromagnetic Variant 

We could consider a very heavy, and hence immobile, nucleus of charge Ze around which 

electrons of mass m and charge — e move. The equations of motion would be identical to 

those in (TTT) but the charges on the r.h.s. replacing the masses and the angular velocity 

becoming now u = — — . 

mR 6 

This idea has been exploited in atomic physics. The most basic example being semi- 
classical treatments of the Helium atom [13]. The idea also extends to muonic atoms H 

C. Time dependence 

Oh et al. |24j . and Heggie et al. [2HJ I2S] mention, in the context of galactic dynamics, a 
time dependent version of ([!]) in which the radius R is allowed to depend on time. In the 
planar case they are 

/.. . /q 2 ^ , o ^ ^ Gm a m b (x b - x a ) 

m a { x a - 2uy a - (3u - 2— )x a + 2uj— y a ) = ) . 

R R tf x a - x 6 \ 6 

b=f=a 

(.. , . R \ Gm a m b (y b - y a ) 
m a {y a + 2ux a - —y) = > . . (63) 

b=f=a 

Generalized Galilei invariance (|7| still holds but (JsJ) becomes 



R R 

x — 2uy — (3u 2 — 2 — )x + 2uj—y = 0, 
R R 



y + 2ux - = . (64) 
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For given R(t) this has a four parameter family of solutions but if one works out the 
Killing vector fields and takes the bracket with time translations d t the algebra will not, 
in the generic case, close on a finite dimensional Lie algebra, even if one adds additional 
generators. The situation is reminiscent to the one considered in Refs. [33H35] . Details will 
be presented elsewhere. 

XI. CONCLUSION 

A remarkable aspect of Hill's equations is that our Eqn. ([8]), simultaneously describes 
time- dependent symmetries ([7]) and the motion of the center of mass. Our solutions ^ 
represent therefore the trajectories both of the symmetry group acting on space-time, and 
of the center-of-mass. 

As long as we consider the 3-body problem it would be physically more important to 
study the relative-motion equation (|6]); the center of mass has little interest, for, say, Lunar 
motions. But Hill's equations also arise when describing an electron beam in a synchrotron; 
guiding center motion is plainly interesting for the latter, as it is in plasma physics, or in 
stellar dynamics [26] . 

As explained in [3H5], the ability to split off the centre of mass motion relies on Galilean 
(in fact Newton-Hooke) symmetry. In our case here rotations are broken, but our time- 
dependent symmetries suffice. 

Moreover, the motion of the center-of-mass can further be decomposed into that of the 
guiding center and relative motion; the generalization of the chiral decomposition [22 | I23), 132] 
is ideally suited for that. It is worth mentioning that, as in [5j 123], our calculations could 
be extended to "exotic" [i.e. non-commutative] particles. 
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